Dielectronic recombination of W 20+ (4d 10 4f 8 ): addressing the half-open f-shell 
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A recent measurement of the dielectronic recombination (DR) of W 20+ [Schippers et al Phys. Rev. 
A 83, 012711 (2011)] found an exceptionally large contribution from near threshold resonances 
(< 1 eV). This still affected the Maxwellian rate coefficient at much higher temperatures. The 
experimental result was found to be a factor 4 or more than that currently in use in the 100 — 
300eV range which is of relevance for modeling magnetic fusion plasmas. We have carried-out DR 
calculations with AUTOSTRUCTURE which include all significant single electron promotions. Our 
intermediate coupling results are more than a factor of 4 larger than our LS-coupling ones at 1 eV 
but still lie a factor 3 below experiment here. If we assume complete (chaotic) mixing of near- 
threshold autoionizing states then our results come into agreement (to within 20%) with experiment 
below < 2 eV. Our total IC Maxwellian rate coefficients are 50-30% smaller than those based-on 
experiment over 100-300 eV. 

PACS numbers: 34.80 Lx, 34.80 Kw 



I. INTRODUCTION 

Tungsten will be a key element [l[ in the ITER mag- 
netic fusion device [2j-|4| currently under construction at 
Cadarachc in France Q- Its ability to withstand high 
power-loads means that it will be the primary facing ma- 
terial within the vacuum vessel. Its high nuclear charge 
means also that it is potentially a serious contaminant in 
the sense of its ability to quench the fusion reaction due 
to radiative power loss. Intensive studies are underway 
at all of the world's major magnetic fusion laboratories to 
understand, predict, and control its behavior. The recent 
ITER-like wall upgrade at JET [f| provides the closest re- 
actor environment short of ITER itself. The 74 ionization 
stages of tungsten may seem daunting from the detailed 
theoretical perspective. Reality is somewhat different. 
Very few ionization stages are observed in practice. Near 
neutrals are seen as they sputter-off surfaces. Then many 
ionization stages burn through quite rapidly before ions 
of much higher charge-state are observed in particular 
localized environments. W 20+ is a significant stage for 
spectral diagnostics and is seen at the null point of the 
separatrix at JET. W 44+ performs a similar function at 
the core and is observed by the JET KX1 spectrometer. 
The ionization stages will change with the much larger 
and hotter ITER device but the principle remains the 
same: very few stages need to be modeled in detail. The 
great bulk of them can be modeled more 'coarsely' as 
supcrstages. Detailed studies are being made of these 
key stages. One of the most basic and important theo- 
retical quantities is the tungsten ionization balance since 
it is the main determinant of the intensity of spectral 
line emission. Electron-collisional equilibrium is a bal- 
ance between electron-impact ionization and dielectronic 
recombination. (All other recombination processes are 
negligible in the magnetic fusion domain.) A sufficiently 
accurate theoretical description of dielectronic recombi- 



nation is key. 

A recent experiment on W 20+ ions by Schippers et al 
Q at the TSR storage ring at Heidelberg measured an 
exceptionally large dielectronic recombination merged- 
beams 'rate coefficient' at a few eV. So-much-so that 
its contribution at the temperatures of significant frac- 
tional abundance for W 20+ (100-300 eV) gave rise to 
a Maxwellian rate coefficient that is a factor 4 or more 
larger than that currently used by the main magnetic fu- 
sion modeling package — the Atomic Data and Analysis 
Structure (ADAS) Q. We seek to resolve this discrep- 
ancy. 

There is little previous detailed work on DR for f- 
shell ions. At one end (4d 10 4/) there are calculations for 
Gd 17+ by Dong et al [l(| utilizing the Flexible Atomic 
Code (FAC) [lj which are relevant for modeling soft 
X-ray lithography. At the other end (4rf 10 4/ 13+ ) there 
are calculations for Au 20+ by Ballance et al [12| with 
AUTOSTRUCTURE []]| . The results of Ballance et al are 
in good agreement with the measurements of Schippers 
et al 0] from 2 meV up to 10 eV. Both of the above ap- 
proaches are standard level-resolved calculations which 
allow for single electron promotions-plus-capture. They 
are largely restricted to at most doubly-excited configu- 
rations and interactions thereof. 

Previous work on open f-shell ions is apparently limited 
to the configuration average approximation, the Burgess 

General Program (bbgp) mm, 

and others of that 
ilk which are the mainstay of modeling codes. The cal- 
culations of Flambaum et al [l6| for the DR of Au 25+ 
(which is isoclcctronic with W 20+ ) can be viewed as a 
form of partitioned configuration average. They utilize 
expressions for the radiative rate and autoionization rate 
which are similar to those of the configuration average. 
The near-threshold autoionization rates are partitioned 
over a Breit-Wigncr distribution which is characterized 
by a single spreading width (l~7| . This compares with 
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our previous configuration average work [l8| which parti- 
tions them over the non-metastable core levels according 
to their statistical weight. Such a partitioning maintains 
the allowed- vs forbidden-channel nature which is char- 
acteristic of DR in simple systems. The justification of 
the Breit-Wigner form and the spreading width follows 
from the complexity of the open f-shell ions in which con- 
figuration mixing tends to a chaotic limit which can be 
described by statistical theory . This leads to a struc- 
tureless continuum for the near-threshold merged-beams 
DR rate coefficient. The Au 25+ results of Flambaum et 
al [l6| were found to be in good ag reement with the 
measurements of Hoffknccht et al [20| below 1 eV. 

The approach to DR reported-on in this paper is a de- 
tailed (level-to-level) one with single electron promotions- 
plus-capture for all significant contributing configura- 
tions. The configuration mixing (and spin-orbit mixing) 
that we allow-for is the same as done in our previous 
work on the open Fe M-shell (3p q ) [2l| and the open 
Sn N-shell (Ad q ) [22[ ions. This allows for the mixing 
of all autoionizing configurations within the (N + 1)- 
electron complex when the promoted and captured elec- 
trons have the same principal quantum number. The 
inequivalent case restricts the mixing to configurations 
of the iV-electron core i.e. with common Rydberg nl 
quantum numbers. This approach gives rise to a near- 
structureless continuum for the near-threshold merged- 
beams DR rate coefficient for these Fe and Sn cases. We 
seek to extend this work to the half-open f-shell and to 
compare the near-threshold merged-beams DR rate coef- 
ficient with the measurements of Schippers et al 0] • We 
seek also to determine the high energy Maxwellian DR 
rate coefficient applicable to the diagnostic modeling of 
W 20+ in magnetic fusion plasmas. 

We provide a description of our background theory in 
Sec. HI1 We describe its application to W 20+ in Sec. UTT1 
We present out results and compare them with those of 
experiment in Scc. lIVl We make some concluding remarks 
in Sec. El 



II. THEORY 

We use the independent processes and isolated reso- 
nance approximations to describe dielectronic recombi- 
nation Il4j. The oft-repeated working equations may be 
found in [2l| along with a more expansive discourse on 
the methodology we employ. We summarize some perti- 
nent points. 



A. Methodology 

We use the computer code autostructure [l3[ to 
calculate all relevant atomic parameters: energy lev- 
els, radiative rates, and autoionization rates. A multi- 
configuration expansion is used in an iS'-coupling or in- 
termediate coupling (IC) representation. The configu- 



ration average (CA) representation [23|, [24| is a simpler 
approach which is very useful for complex heavy species 
since it provides a rapid overview of the problem. 

All of our previous works with the CA approxima- 
tion utilized the dracula code [H, [Hi- One im- 
mediately sees the effect of level-resolution and fine- 
structure mixing when comparing L^-coupling with IC 
results obtained with autostructure. Comparison 
of LS'-couping with CA results has been clouded by 
the fact that dracula is based upon the Cowan struc- 
ture code [2(|. The Cowan code utilizes kappa-averaged 
Hartree-Fock radial orbitals [27j . The differences here 
with autostructure can be minimized by its use of 
kappa-averaged orbitals computed in self-consistent con- 
figuration average model potentials |2f|. The Cowan 
code generally also scales various operator interactions 1 . 
This facility is not readily or generally available in 
autostructure. We have implemented the CA angu- 
lar momentum representation within AUTOSTRUCTURE. 
This eliminates any uncertainty in seeing the pure effect 
of moving to a configuration-mixed term-resolved repre- 
sentation. We have carried-out detailed CA comparisons 
between autostructure and DRACULA in the course of 
the present work so as to verify the integrity of the new 
development. 



B. Computation 

The (near-) half-open f-shell problem is a daunting one. 
If we view it simply in terms of binomial coefficients for 
the number of states present in a configuration then mov- 
ing from the 4d 10 4/ 13 ground configuration considered 
in [l2j to 4cf 10 4/ 8 increases the number of states by a 
factor (429/2). Memory requirements (CPU and disk) 
scale as (429/2) 2 and the time requirement as (429/2) 3 . 
This is all relative. Absolute numbers are much larger 
once we start promoting electrons from the 4d— and 4/— 
subshells of the ground configuration. (The only bonus 
of the binomial effect is that 4/ 7 is only marginally worse 
than 4/ 8 since the number of states increases only by a 
factor 8/7.) 

The published autostructure code [H[ was used 
for our recent work on the tin half-open d-shell [22J and 
Au 20+ 4/ 13 [12]. It does not scale to the half-open f-shell. 
Substantial development has been necessary. A detailed 
exposition is more suited to a computer physics journal 
and so we give only a flavor here. 

Most angular momentum packages used by atomic 
structure and collision codes are based upon Racah al- 
gebra. This is a hierarchical coupling scheme. A com- 
plication for the open f-shell is the need to introduce a 
new quantum number: seniority, autostructure em- 
ploys the non-hierarchical Slater-state approach to angu- 



It does not do so under CA operation. 
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lar momentum coupling that was advocated by Condon 
and Shortley [28|. It has no concept of parentage 2 . All 
interactions are expanded and determined initially in an 
uncoupled (Slater-state) representation. A transforma- 
tion is then made to an LS and/or LSJ representation 
through the use of vector-coupling coefficients (VCCs). 
The half-open f-shell requires billions of coefficients and 
so tens of Gbytcs of RAM per processor. They are not 
all required at the same time. It has been possible to 
implement an archival on disk in such a way that does 
not swamp the calculation with I/O. 

Another issue concerns the large number of radiative 
and autoionization rates that arise in a level-resolved cal- 
culation. The approach to-date has been to archive them 
all to disk for subsequent processing in a variety of ways: 
to compare with experiment, to generate total rate coef- 
ficients for astrophysical modeling, or to process as final- 
state resolved data for collisional-radiative modeling for 
magnetic fusion. Such general flexibility comes at the 
cost of many Tbytes of disk space and corresponding I/O 
time. If we sacrifice some degree of generality then we can 
carry-out bundling over quantum numbers and summa- 
tion of partial widths on-the-fly as the atomic data is gen- 
erated. This reduces the data files to a manageable size. 
The user choice of bundling and/or summation should be 
guided by the exact same implementation made within 
the collisional-radiative modeling approach so as to ren- 
der it tractable for heavy species. It is important to note 
that this introduces no additional approximation for our 
description of the experimental DR cross section nor the 
total Maxwellian rate coefficients presented later. 

We do not dwell on various RAM and CPU time issues 
that arose as well on scaling from the open d- to f-shell. It 
is sufficient to note that the calculations reported-on be- 
low took about two weeks on a modest cluster (the prob- 
lem does not engender large-scale parallelization) with 
4Gb RAM and 250Gb of scratch disk per processor. 



C. Mixing 

We discuss the role of (configuration) mixing of au- 
toionizing states on DR in complex heavy species. This 
is important since we include only a limited amount (see 
the next section). Consider a model problem in which 
we include only autoionizing configurations which re- 
sult from single electron promotions-plus-capture from 
the the ground configuration. We assume that the au- 
toionization rates (Ai) and radiative rates (R) satisfy 
Ai » R. Then the DR cross section is basically pro- 
portional to R for a fixed symmetry. (We consider only 
autoionizations which are inverses of the dielectronic cap- 
ture in this model i.e. near-threshold.) Now consider the 



2 It is both possible and advantageous to introduce some parentage 
but it is not required. 



addition of autoionizing configurations which result from 
double electron promotions-plus-capture from the ground 
configuration. This set is only populated by mixing with 
the first via a unitary transformation. We denote the 
second set of autoionization rates by Ai and assume sim- 
ilar radiative rates (R). We assume A2 << R. If mix- 
ing between the two sets is strong enough and there are 
enough states then the autoionization rates for the first 
set (which we label now as A12) are depleted to the ex- 
tent that A12 << R. The DR cross section is now pro- 
portional to A12 + A2 = A\. The total is enhanced by 
a factor A\/R. (Mixing only takes place between states 
of the same statistical weight and so the sum over au- 
toionizing states is assumed to be implicit.) We con- 
sider the effect of further mixing. Now add configurations 
which result from triple electron promotions-plus-capturc 
from the ground configuration. Denote the autoioniza- 
tion rates for the three sets as A123, A23, and A3. As- 
sume all Ai << R still. The total DR cross section is un- 
changed since it is proportional to A123 + ^23 + A 3 = A\ 
still. Such mixing is merely redistributive and so can be 
neglected with respect to total DR. 

This is the nature of the near-threshold DR pro blem 
in complex ions described by Flambaum et al 11611 (for 
Au 25+ ) following a preliminary earlier study [171 ]. The 
chaotic fully-mixed nature of Au 24+ was verified by Grib- 
akin and Sahoo [29[ who performed large-scale calcula- 
tions of eigen-encrgies and eigen- vectors. Our detailed 
description of the autoionizing spectrum is necessarily 
incomplete. The above discussion is intended to shed 
light on how far we need to go. The model problem was 
merely illustrative in using configuration mixing via dif- 
ferent sets of promotions. We are largely restricted to one 
electron promotions-plus-capture since we need to com- 
pute autoionization rates and radiative rates which are 
applicable over a much wider range of energies. But we 
do have differing representations: CA, LS'-coupling and 
IC. The question then is one of the degree to which the 
IC representation is incomplete with regards to enhance- 
ment or whether it has moved to a redistributive regime 
so that our DR cross sections will have converged largely 
with respect to the total. The structured behavior of 
our low energy cross section, its comparison with experi- 
ment, and with statistical theory will enable us to judge 
the degree to which this is so. (It is straightforward to 
apply the statistical approach to our data — we simply 
partition our autoionization rates using the Breit-Wigner 
distribution — but note that this procedure is only valid 
near threshold.) 



III. APPLICATION 

The ground configuration of W 20+ is [Kr]4d 10 4/ 8 . The 
ground term is 7 F. The ground level is J = 6. These 
denote the ground state of the CA, LS'-coupling, and 
IC representations. We describe the DR reactions that 
we take account-of by their configuration representation. 
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This consists of TV-electron configurations to which a con- 
tinuum electron and a Rydberg electron are each cou- 
pled. This describes dielectronic capture and/or au- 
toionization. The latter describes radiative transitions 
within the core also. Additional (N + l)-electron 'cor- 
relation' configurations are added to describe Rydberg 
electron radiation into the core. Rydberg-Rydberg ra- 
diation n — > n' > A is described hydrogcnically. We 
breakdown the problem by target subshell promotion. 



A. An = 

We allow for Ad —> 4/ and 4/ — s> 4/ promotions in 
both L5-coupling and IC. The latter promotion does not 
contribute to CA DR since it corresponds to an clastic 
transition. The TV-electron configurations are 4c? 10 4/ 8 
and 4<i 9 4/ 9 . The (N + l)-electron configurations are 
4d 10 4/ 9 and 4<i 9 4/ 10 . We note that some terms/levels 
of Ad 10 A f 7 51 lie below 4<i 9 4/ 9 . The former could pro- 
vide an alternative autoionization channel for the lat- 
ter. Such a transition is forbidden directly since it is 
described by a two-body operator. (A continuum elec- 
tron is still to be coupled to the former and a Rydberg 
to the latter.) It could take place via mixing if we were 
to include additional configurations. We do not. We 
consider also the Ap — > 4/ promotion in CA only. Its 
contribution is expected to be small. We normally use 
the CA calculation to determine the range of Rydberg- 
nl required to converge the total DR and then use these 
values subsequently for the more demanding LS'-coupling 
and IC calculations. We find that the contribution to the 
Maxwcllian rate coefficient from the 4c? — > Af promotion 
is converged to about 3% by n = 100 and I = 7 at all 
energies. It is not possible to do so for the 4/ — > 4/ pro- 
motion since the CA result is zero and so we utilize the 
LS'-coupling calculation here to delimit the IC one. We 
find that the contribution from the 4/ — > 4/ promotion 
is converged to about 3% by n = 100 and I = 6 at all 
energies. 



B. An=l 



We consider the Ad and Af promotions separately. 



1. Af ->■ 5/ 

The iV-electron configurations are 4/ 8 and Af 7 5l (I = 
— 4). It is necessary to omit the n = 5 continuum so as 
to keep the problem tractable in the L5-coupling and IC 
calculations 3 . We carried-out CA calculations both with 
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FIG. 1: (Color online) W 20+ CA Maxwellian DR rate coef- 
ficient contributions for various promotions: total (solid red 
curve), Ad — > Af (long-dashed green curve), Af — > 51 (short- 
dashed blue curve), Ad — > 51 (dotted magenta curve), Ap — > Af 
(dot-dashed cyan curve) , and Ad + Af — > 61 (double-dashed 
orange curve). 



and without the n = 5 continuum to aid our analysis of 
the uncertainty (overestimate) in our LS'-coupling and 
IC results. There is none at all below ~ 20 eV since they 
are all energetically closed. 

The (N + l)-electron configurations are Af s 5l (I = 

— 4). Some of the Af 7 5l5l' configurations arc (partially) 
bound. We treat such n — > n' = 5 radiation hydrogeni- 
cally (approximately) in the LS'-coupling and IC calcu- 
lations. They are either strictly bound or autoionizing in 
the CA approximation. Their contribution is small. 

The CA results for this promotion are converged to 
about 5% at 200eV and 10% at 1000 eV on summing to 

1 = 5. The sum over n is converged to better than 2% 
by n = 100. We add this small 'top-up' in I (and n since 
we are doing so) from our CA results to the coupling 
and IC ones. 



3 The number of VCCs that need to be internally buffered becomes 
too large. A smaller buffer could be implemented but this would 



likely increase I/O time substantially. The absolute number of 
VCCs required in LS-coupling is only typically a factor of 2 
smaller than for IC. The demands of the IC calculation arise 
from the fact that far more of the states that they represent 
interact. 
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2. 4d -> 5/ 

The TV-electron configurations are 4d 10 4/ 8 , 4d 9 4/ 9 , 
4ci 9 4/ 8 5Z (1 = 0-4) and id 10 4 f 7 51 (1 = 0-4). We omit 
the n — 5 continuum again. The 4d 9 4/ 8 is rather de- 
manding when coupled to 5lnV for I + V > 5. It has a 
factor 70/8 more states than the corresponding 4d 10 4f 7 . 
We need to consider it further. We write the dielectronic 
capture reaction in a somewhat unusual form: 



10000 



4f( 7 F)4d 1 



4f( 7 F)4d 9 5lnl' . 



This illustrates the role of the 4/ 8 7 F ground term as a 
spectator. It cannot change simultaneously with the two- 
body dielectronic capture. It can change (shake-up) via 
mixing in the autoionizing states. We omit such mixing. 
We do the same for the reverse radiative stabilization to 
4/ 8 ( 7 F)4d nl. This renders a tractable but reasonable 
description of the 4c? — > 51 promotions. We recall that 
there is no configuration mixing whatsoever present in 
the CA approximation. We recall also the lack of sensitiv- 
ity to configuration mixing that we found for total DR in 
the tin 4p — 4d transition array despite the demonstrable 
extensive configuration mixing [22|. Such an argument 
may not be valid near threshold — see the discussion in 
Sec. Ill CI We emphasize that we place no such restric- 
tions (on mixing) when the 4/ is active such as for the 
4d 9 4f 9 configuration. We have implemented the general 
user specification of term restrictions for spectator sub- 
shells within AUTOSTRUCTURE. These are common to all 
configurations which contain the specified subshcll(s). 

The (N + l)-electron configurations are 4<i 10 4/ 8 5Z (I = 
0-4) and 4d 9 4f 9 5l (1 = 0-4). A few of the 4d 9 4f s 5l5l' 
configurations are (partially) bound and we treat them 
as for 4/ ->• 51. 

The CA results for this promotion are converged to 
1% at 1000 eV on summing to I = 4. The sum over n is 
converged to much better than 1% by n = 100. 



C. A?i = 2 

We consider 4d + 4f — » 61 promotions within the CA 
approximation only. The contribution is expected to be 
small. 
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FIG. 2: (Color online) W 20+ Maxwellian DR rate coefficient 
contributions for An = promotions: IC 4d — s> 4/ (solid 
red curve), IC 4/ — >• 4/ (long-dashed green curve), LS 4d — ► 
4/ (short-dashed blue curve), LS 4/ — > 4/ (dotted magenta 
curve), and CA 4d — > 4/ (dot-dashed cyan curve). 



We see that we do not need to consider 4p — > 4/ and 
An = 2 promotions any further. 



A. An = 

We present and compare our LS-coupling and IC re- 
sults for the 4/ — > 4/ and 4c! — s> 4/ promotions in Fig. [2] 
We note the close agreement between them and the CA 
results for the 4d — > 4f promotion above ~ 100 eV. There 
is about a factor 7 difference between the L5-coupling 
and IC results down at 1 eV. The contribution from the 
4/ — > 4/ is no more than about 10% of the 4d -> 4/ 
above 1 eV. 



IV. RESULTS 



We show an overview of the different CA contributions 
to the total DR Maxwellian rate coefficient in Fig. [TJ The 
energy range of interest for an electron collisional plasma 
is 100-300 eV. This is where W 20+ has its maximal frac- 
tional abundance (> 0.01) in a magnetic fusion plasma. 

We remark that the contribution from 4d — > 4/ pro- 
motions is comparable with the An = 1 above ~ 100 eV. 
This is in contrast to case of the almost full 4f-subshell 
case of Au 20+ [l2j. It also dominates at a few eV but 
this behavior can be expected to be more ion-dependent. 



B. An = 1 

1. 4/ -»■ 51 

We present and compare our L5-coupling and IC re- 
sults for the 4/ — > 51 promotions in Fig. [3] We separate- 
out the contributions from capture to n = 5 and n > 5. 
The sum of the two in LS'-coupling and IC agree to wit hin 
10% by 100 eV. We show CA results both with and with- 
out autoionization to the n = 5 continuum. The rather 
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FIG. 3: (Color online) W 20+ Maxwellian DR rate coefficient 
contributions for 4/ — > 51 (I — — 4) promotions (capture to 
n = 5 and n > 5 are shown separately) : IC (solid red curves) , 
LS (long-dashed green curves), and CA, n > 5 both with and 
without the n = 5 continuum (short-dashed blue curves). 
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FIG. 4: (Color online) W 20+ Maxwellian DR rate coefficient 
contributions for 4rf — > 51 (I = — 4) promotions (capture to 
n — 5 and n > 5 are shown separately): IC (solid red curves), 
LS (long-dashed green curves), and CA, n > 5 both with and 
without the n = 5 continuum (short-dashed blue curves). 



pronounced high energy peak is reduced by a factor of 
2 at 160 eV. These are the first autoionizations into ex- 
cited states pathways in the CA coupling scheme. The 
LS'-coupling and IC are already suppressed by autoion- 
ization into (the continuum of) a multitude of excited 
states within the ground configuration. This is reflected 
in their less pronounced high energy peaks. We would 
not expect the (n > 5) LS-coupling and IC results to be 
suppressed further by more than ^20% below 300 eV. 



2. 4d ->■ 5/ 

We present and compare our LS-coupling and IC re- 
sults for the Ad — > 51 promotions in Fig. [4] We see that 
the relative contributions from capture to n = 5 and 
n > 5 are reversed compared to the 4/ — > 5/ promo- 
tions. The n > 5 contribution does not exceed that of 
the n = 5 until high energy. This is due to autoionization 
suppression via the 4/ — > Ad inner-shell transition. The 
CA results for n > 5 are suppressed by a factor of two at 
160 eV but the sum including n = 5 by about one third. 
The (n > 5) LS'-coupling and IC results are likely to be 
suppressed by a larger relative factor than for the case of 
4/ — > 51 promotions but the overall sum including n = 5 
dilutes the factor and ~ 20% appears to be a reasonable 



estimate here. 



C. Totals (Merged-beams) 

We convoluted our DR cross sections with the electron 
velocity distribution applicable for the merged electron- 
ion beams in the TSR cooler Q. We compare our re- 
sultant total DR rate coefficients with those measured 
on the TSR storage ring Q- We focus first on the near 
threshold region [0-10] eV for which Schippers et al re- 
ported the largest measured DR rate coefficient to-datc. 
We see in Fig.[5]that our IC results are a factor 2-5 larger 
than our LS'-coupling ones but they are still typically a 
factor of 3 smaller than experiment. Our IC results are 
dominated by the id — > 4/ promotion (80%). The den- 
sity of resonances is such that there is little resultant 
structure in the total. We have a near quasi-continuum 
of resonances. It does not matter then just where the 
ionization limit lies. 

The remaining factor of ~ 3 is likely due to the in- 
complete mixing within our IC configuration expansion. 
We show results (Fig. [5]) where we have partitioned our 
autoionization rates using the Breit-Wigner distribution 
with a spreading width of 10 eV [17[. The results are not 
particular sensitive to this width. Our CA, LS'-coupling 
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FIG. 5: (Color online) W 20+ merged-beams DR rate coeffi- 
cients: experiment [?J (solid black curve), partitioned total 
(dot-dashed cyan curve), IC total (solid red curve), LS to- 
tal (long-dashed green curve), and IC 4d — > 4/ only (short- 
dashed blue curve). 



and IC results are barely distinguishable on this scale 
and so we show a single curve. We obtain agreement 
with experiment to within 20% over 0.003 - 2 eV. Sim- 
ilar findings were obtained by Flambaum et al for 
Au 25+ . (The measured cross section increases at energies 
below 0.003 eV due to an artifact of the merged-beams 
technique f3Cj.) 

The first excited level of the ground term opens-up just 
below 2 eV. The experimental cross section falls away 
progressively thereafter — see Fig. [5J This fall-off coin- 
cides with an increasing number of alternative autoion- 
ization channels opening-up. If autoionization into ex- 
cited levels is fully redistributed as per the ground then 
the total width is unaffected. The partitioned autoion- 
ization widths to the ground level are orders of magni- 
tude smaller than the radiative widths. This is why the 
partitioned results are largely unchanged — even when 
summing over autoionization to hundreds of excited lev- 
els. In nature it would appear that typical autoionizing 
widths to the ground level are no more than a factor 
~ 10 smaller than the radiative widths. Summing over 
a relatively small number of excited continua then pro- 
duces a total autoionization width comparable with and 
then exceeding the radiative width. We note that our 
level-resolved autoionization widths are only typically at 
most a factor 5-10 larger than the radiative widths. Ex- 




25 50 75 100 125 
Energy (eV) 

FIG. 6: (Color online) W 20+ merged-beams DR rate coeffi- 
cients: experiment [?J (solid black curve), partitioned total 
(dot-dashed cyan curve), IC total (solid red curve). 



pcrimcnt falls between the two theoretical 'limits'. The 
partitioned results are clearly inapplicable here. 

All of our results are an upper limit because we assume 
100% of the W 20+ initial ion population to be residing 
in the ground state. Schippers et al Q identify several 
possible metastable levels that could remain populated 
during the lifetime of the measurement but have no esti- 
mate of their population. Any combination of mctasta- 
bles that we take reduces the total. This is because DR 
from excited states is suppressed by autoionization into 
the continuum attached to lower levels. Only the ground 
level is immune from such. The agreement between our 
Breit-Wigner partitioned results and experiment below 
10 eV would indicate that the metastable presence is not 
significant. 



D. Totals (Maxwellian) 

We turn now to the corresponding total Maxwellian 
DR rate coefficients. We compare results in Figs. [7] and 
[51 The two CA results (Fig. [7J illustrate the effect of 
omission of n = 5 continuum suppression on the total. 
It is 25% at 160 eV. The lower CA, LS'-coupling, and IC 
results are all in close agreement above ~100 eV. The 
LS'-coupling and IC results are an upper limit not just 
because of the omission of the n = 5 continuum but also 
because they assume 100% of the W 20+ initial ion popu- 
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FIG. 7: (Color online) W 20+ total Maxwellian DR rate coef- 
ficients: IC (solid red curve), LS (long-dashed green curve), 
and CA with-and-without n — 5 continuum (short-dashed 
blue curves). The fractional abundance of W 2 
netic fusion plasma is shown also (solid black curve). 
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lation to be residing in the ground state. 

The experiment by Schippers et al Q only detects 
resonances which occur below 140 eV. We show a second 
IC result (Fig.© which imposes such a cut-off. This cut- 
off result lies just over 50% below experiment at 160 eV. 
If we add a theoretical top-up (for the resonances above 
140 eV) to the experimental result then this difference is 
reduced to about 40% (at 160 eV) and is 50-30% over 
100-300 cV. The top-up is between 20-60% of the ex- 
perimental result alone over 100-300 cV. The topped-up 
experimental result is the total (zero-density) DR rate 
coefficient that we recommend for use in modeling be- 
cause of the remaining difference between theory and 
experiment. This difference is attributable to the dif- 
ference in the contribution from low energy resonances. 
The An = resonances as a whole contribute about one 
third of the IC total at 160 cV. The factor 3 larger experi- 
mental contribution from such resonances at low energies 
means that they contribute significantly more here. 

The final results we show in Fig. [8] are those from 
ADAS [3ll |. These were determined using the Burgess 
General Formula (GF) [H — this is ADAS Case A @ 
which extends the GF validity to finite density by the 
use of a global suppression factor (33|. Those shown 
here were determined at close to zero electron density 
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FIG. 8: (Color online) W 20+ total Maxwellian DR rate co- 
efficients: IC all resonances and to 140 eV only (solid red 
curves), experiment [?J to 140 eV and with theory top-up 
for resonances above 140eV (long-dashed green curves), and 
ADAS [U (short-dashed blue curve). The fractional abun- 
dance of W 20+ in a magnetic fusion plasma is shown also 
(solid black curve). 



We see that the ADAS [3i( results lie between a factor 
10 and 4.5 below our recommended ones over 100-300 cV. 
The ADAS results are for dielectronic plus radiative re- 
combination. The radiative contribution dominates in 
these results below 20 eV because the Case A Burgess 
GF cannot describe the effect of low-lying resonances. 
Such low-lying resonances can be described by the the 
BBGP approach [3, [H] — this is ADAS Case B which 
resolves the final recombincd state and so is amenable to 
the collisional-radiative modeling of density effects [34| • 

The fractional abundance curve we show has been de- 
termined using the ADAS Case A data at an electron 
density of 10 13 cm~ 3 which is typical of that relevant to 
magnetic fusion edge plasmas [l[ . It differs slightly from 
the one shown by Schippers et al Q which is due to 
Putterich [35| and at 10 14 cm~ 3 . It is appropriate to use 
a finite-density abundance to indicate the plasma rele- 
vant temperatures on which to focus our comparisons of 
zero-density rate coefficients because the temperature of 
peak abundance is density sensitive. Use of rate coeffi- 
cients which are up to a factor 10 larger though is likely to 
move the peak abundance to higher temperature. (Simi- 
lar increases in the DR rates to used can be expected for 
adjacent ionization stages.) 
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TABLE I: Recommended total (zero-density) dielec- 
tronic recombination rate coefficient fitting coefficients 
d (cm 3 s _1 [eV] 3 ' 2 ) and Ei(eV) for the initial ground level of 
W 20+ . 



i Ci Ei 

1 4.025(-7)° 1.093(+0) 

2 7.697(-7) 9.153(+0) 

3 1.065(-6) 3.425(+l) 

4 1.487(-6) 1.205(+2) 

5 2.177(-6) 2.384(+2) 



a (m) denotes x!0 : 



In Tabic U we present such fitting coefficients for the 
recommended (experiment topped-up by theory) total 
zero-density dielectronic recombination rate coefficient 
for the initial ground level of W 20+ . The fit is accurate 
to better (often much better) than 1% over 1-1000 eV. 

The total dielectronic recombination rate coefficient 
can be taken to be the total recombination rate coeffi- 
cient. The contribution from radiative recombination is 
negligible over the given temperature range as is that 
from three-body recombination at the densities of inter- 
est to magnetic fusion plasmas. 



1. Density effects 

A rigorous treatment of density effects on dielectronic 
recombination rate coefficients and their consequential 
effect on the ionization balance of W 20+ and adjacent 
ionization stages is beyond the scope of the present work. 
We can make some observations though. The ADAS rate 
coefficient is reduced by a factor of 2 at an electron den- 
sity of 10 13 cm~ 3 (not shown) compared to zero-density. 
This is due to the stepwise ionization of high-n (> 10) 
Rydberg states following recombination. The new recom- 
mended total DR rate coefficients contain a large contri- 
bution from low-energy resonances of low-n (< 10). If we 
assume that the high-n contributions to both are broadly 
similar then we can expect maybe a 10-20% reduction 
in the new recommended values over 100-300 eV. Simi- 
lar (reduced) effects can be expected for adjacent ioniza- 
tion stages. This means that the corresponding fractional 
abundances are likewise less sensitive to the electron den- 
sity than indicated by the current ADAS data. A revi- 
sion of the density dependent ionization balance of f-shcll 
tungsten ions is clearly needed. 

2. Fitting coefficients 



V. CONCLUSION 

We have calculated IC DR rate coefficients for W 20+ 
which include all significant one-electron promotions- 
plus-capture. A factor 3 difference with experiment re- 
mains at low energies. We have demonstrated that this 
can be removed if wc assume complete chaotic mixing of 
multiply-excited near-threshold configurations. A similar 
finding was obtained by Flambaum et al [l6| for Au 25+ . 
The difference between theory and (topped-up) experi- 
ment at energies relevant to magnetic fusion modeling for 
ITER is somewhat less viz. between a factor two and 1.5 
over 100-300 cV. The DR data used by ADAS for such 
modeling needs to be updated since the current Burgess 
GF Case A results lie between a factor 10 and 4.5 be- 
low our new recommended values over 100-300 cV. Our 
CA, iS-coupling and IC results are all in close accord 
(20%) above lOOeV which suggests that DR rate coeffi- 
cients for complex W ions can be determined readily to 
within a factor of two for modeling purposes. Similar be- 
havior can be expected for related complex ions of other 
heavy elements. Determination of such DR rate coef- 
ficients which are accurate to the ~ 20% level remains 
problematic though. 



It is convenient often for simple modeling purposes to 
fit the total Maxwellian dielectronic recombination rate 
coefficient (a) to the following functional form: 

a(T) = T- 3 / 2 5>cxp(^i) 

where the Ei are in the units of temperature T (e.g. eV) 
and the units of Ci are then cm 3 s _1 [eV] 3 / 2 . 
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